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Abstract 

We report on a comprehensive study of the phase structure of cold, dilute nuclear matter fea- 
turing a 3 S±- 3 Di condensate at non-zero isospin asymmetry, within wide ranges of temperatures 
and densities. We find a rich phase diagram comprising three superfluid phases, namely a LOFF 
phase, the ordinary BCS phase, and a heterogeneous, phase-separated BCS phase, with associated 
crossovers from the latter two phases to a homogeneous or phase-separated Bose-Einstein conden- 
sate of deuterons. The phase diagram contains two tri-critical points (one a Lifshitz point), which 
may degenerate into a single tetra-critical point for some degree of isospin asymmetry. 
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Introduction. Fermionic Bardeen-Cooper-Schrieffer (BCS) superfluids, which form loosely 
bound Cooper pairs at weak coupling, undergo a transition to the Bose-Einstein condensate 
(BEC) state of tightly bound bosonic dimers, once the pairing strength increases jl, 2|. This 
behavior has been confirmed in experiments on cold atomic gases, where the interactions can 
be manipulated via the Feshbach mechanism In isospin-symmetric nuclear matter, 

the transition from the BCS to the BEC state of the 3 Si- 3 -Di condensate may occur upon 
dilution of the system, in which case the asymptotic state is a Bose-Einstein condensate of 
deuterons Js - 12] . 



Isospin asymmetry, induced by weak interactions in stellar environments and expected 
in exotic nuclei, disrupts isoscalar neutron-proton (np) pairing, since the mismatch in the 



Fermi surfaces of protons and neutrons suppresses the pairing correlations [13|. The standard 
Nozieres-Schmitt-Rink theory j^j of the BCS-BEC crossover must also be modified, such that 
the low-density asymptotic state becomes a gaseous mixture of neutrons and deuterons f^ |. 

Two relevant energy scales for the problem domain under study are provided by the 
shift 5fi = (fi n — /i p )/2 in the chemical potentials \x n and \i v of neurons and protons from 
their common value \x§ and the pairing gap Ao in the 3 Si- 3 Di channel at 5fi = 0. With 
increasing isospin symmetry, i.e., as 5fi increases from zero to values of order A , a se- 
quence of unconventional phases may emerge. One of these is a neutron-proton condensate 



whose Cooper pairs have non-zero center-of-mass (CM) momentum 10|, ]l5|, ]l6| ; this phase is 
the analogue of the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase in electronic supercon- 



ductors 
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181 ] . Another possibility is phase separation into superconducting and normal 
components, proposed in the context of cold atomic gases [li|. At large isospin asymmetry, 
where 3 Si- 3 Di pairing is strongly suppressed, a BCS-BEC crossover may also occur in the 
isotriplet l So pairing channel, notably in neutron-rich systems and halo nuclei 2CH27|. 

Our main objective is to combine the ideas of unconventional 3 S\- 3 Di pairing and the 
BCS-BEC crossover in a model of isospin-asymmetric nuclear matter and construct a phase 
diagram for superfluid nuclear matter over wide ranges of density, temperature, and isospin 
asymmetry, while also including non-BCS pairings. By doing so, we advance the computa- 
tional treatment of dilute hadronic matter along several lines, (i) The BCS-BEC crossover 
in isospin-asymmetric systems, studied previously in Ref. [l^], is extended to include, for 
the first time, a phase with broken spatial symmetry and a spatially symmetric but hetero- 



geneous phase, (ii) We extend the earlier studies 
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Q of the nuclear LOFF phase to 
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the low density regime and show that this phase is succeeded by a less dense heterogeneous 
phase before a transition to the BEC regime occurs, (iii) We provide the first treatment of a 
heterogeneous (phase-separated) neutron-proton condensate in the context of 3 Si- 3 -Di-paired 
nuclear matter. 

Theory. The Green's function of the superfluid, written in the Nambu-Gorkov basis, is 
given by 

where Gf 2 = G^Jxi,x 2 ), etc., x = (t,r) denotes the continuous temporal-spatial variables, 
and Greek indices label discrete spin and isospin variables. Each operator in Eq. ([1]) can 
be viewed as a bi-spinor, i.e., ip a = (ip n ^,ip n i,ip p ^,ip P i) T , where the internal variables f, 4- 
label a particle's spin, and n,p its isospin. The matrix propagator obeys the familiar Dyson 
equation, which has the formal solution 



'0,13 



L3 " Sl 3 ) %2 = S 12 , (2) 



in terms of the matrix self-energy S, where summation/integration over repeated indices is 
implicit. Eq. (J2} is advantageously transformed into momentum space, where it becomes 
an algebraic equation. For our purposes, translational invariance cannot be assumed, so 
we proceed by defining relative and center-of-mass (CM) coordinates f = (x\ — x 2 ) and 
R = (xi + x 2 )/2 and Fourier transforming with respect to the relative four-coordinate and 
CM three-coordinate R. The associated relative momentum is denoted below by k = (ik u , k) 
and the three-momentum of the CM is denoted by Q. The zero component of the vector 
k takes discrete values k v = (2u + l)irT, where v e Z and T is the temperature. Thus the 
Fourier image of Eq. (J2J) is written as 

[%(k, Q)- 1 - E{k, Qj] &(h, Q) = l 8x8 . (3) 

The normal propagators for the particles and holes are diagonal in both spaces, i.e., 
(G + ,G~) oc Saa'] hence the off-diagonal elements of ^ _1 are zero. Writing out the non- 
vanishing components in the Nambu-Gorkov space explicitly, we obtain [Wo(ik u , k, Q^^n = 
-[%(-ik u , k, -Q) _1 ] 2 2 = Gq 1 ^, k, Q), where 

G {k)~ l = diag(z&^ - e+ t , %k v - e+ 4 , ik v - ej, ik v - e+ ; ) (4) 
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with e% = e± =E s -5fi±E A and e± = e± = E s + 8fi±E A . Here E s = (Q 2 /4 + k 2 ) /2m* - 
/i and E A = k ■ Q/2m*, with Jx = (/x n + fi p )/2. The effective mass m* is defined in 
the usual fashion in terms of the normal self-energy, bare mass m, and Fermi momentum 
Pf, i.e., m/m* = [1 — (m/p)d p 'E 1 i\ p=PF ], the small mismatch between neutron and proton 
effective masses being neglected. The quasiparticle spectra in Eq. (j3J) are written in a general 
reference frame moving with the CM momentum Q with respect to a laboratory frame at 
rest. The spectrum of quasiparticles is seen to be two-fold degenerate, i.e., the SU(4) Wigner 
symmetry of the unpaired state is broken down to spin SU(2). Note that this symmetry 
is always approximate, since the phase shifts in the isoscalar and isotriplet 5-waves differ, 
such that isosinglet pairing is stronger than isotriplet pairing in bulk nuclear matter. 

The nucleon-nucleon scattering data show that the dominant attractive interaction in 
low-density nuclear matter is the 3 Si- 3 .Di-partial wave, which leads to isoscalar (neutron- 
proton) spin-triplet pairing. Accordingly, the anomalous propagators have the property 
(-^12,-^12) oc (—i<Ty) <8> t x , where <Tj and T{ are Pauli matrices in spin and isospin spaces. 
This implies that in the quasiparticle approximation, the self-energy H has only off-diagonal 
elements in the Nambu-Gorkov space. Specifically, S12 = = iA a p, with A 14 = A23 = 
— A 32 = — A41 = iA, where A is the (scalar) pairing gap in the 3 Si- 3 -Di channel. Substituting 
Eq. (J4]) into Eq. ([21), we obtain a set of algebraic equations whose solutions provide the 
"normal" and anomalous Green's functions 

n/p " (iku ~ E+ /± ){ik v + E-^Y 1 ) 

np (iK-E+wK + E-y [) 

F± = - (7) 

the four branches of the quasiparticle spectra being given by 



E a r = y/E% + A 2 + r5fi + aE A , (8) 

in which a,r G {+, — }. Analytic continuation of these Green's functions via ih v — > k + i0 + 
yields their retarded counterparts. The densities of neutrons and protons in any of the 
superfluid states are obtained through 

d 4 k 



Pn/p 



—ImiG+ /p (k,Q)-G- /p (k,Q)]f(u J ), 

(9) 



where k = (ho, k) and f(x) = l/[exp {x/T) + V\. In mean-field approximation, the anomalous 
self-energy (pairing-gap) is determined by 

A(k,Q) = 2* J0pV(k,hT) 

x Im[F+(k, Q) - F~ n (k, Q)]f(u), (10) 

where V(k,k') is the neutron-proton interaction potential. Performing a partial-wave ex- 
pansion in Eq. ({10"]) as well as the integration over ho, we find 

73 u 



n. v * ' 



x __^E^^=[l-2/(^)], (11) 
2v/^(^) 2 + A,(^Q) 



where Vi^(k,k') is the interaction in the 3 Si- 3 .Di partial wave. The magnitude Q of the 
CM momentum in Eqs. © and (ITT]) is a parameter to be determined by minimizing the 
free energy of the system. For the homogeneous (but possibly translationally non-invariant) 
cases it suffices to find the minimum of the free energy of the superfluid (S) or unpaired (N) 
phase, 

Fs = E s — TS S , F N = E N — TS N , (12) 

where E is the internal energy (statistical average of the system Hamiltonian) and S denotes 
the entropy. Stability of the superfluid phase requires Fs < F^. Three possibilities exist 
for the homogeneous phases: (i) Q = 0, A ^ (BCS phase), (ii) Q ^ 0, A ^ (LOFF 
phase), and A = 0, Q = (unpaired phase). The free energy of the heterogeneous phase 
(phase-separation case) is constructed as a linear combination of the superfluid and unpaired 
energies, 



&(x, a) = (1 - x)F s (a = 0) + xF N (a ^ 0), (13) 

where x is the filling fraction of the unpaired component and a = (p n — p p )/(p n + p p ) is 
the density asymmetry. In the superfluid phase (S) one has p^ = p p = p^ s \ whereas 
in the unpaired phase (N) they are rescaled to new values p^/p- Thus, the net densities of 
neutrons/protons per unit volume are given by p n / p = (1 — x)p^ + xp^j\ 

Results. Eqs. ([9]) and (fTTl) were solved self-consistently for a pairing interaction given by 
the bare nucleon-nucleon interaction in the 3 S\- 3 Di partial wave, based on the (phase-shift 
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equivalent) Paris potential 28j. The nuclear mean-field was modeled by a Skyrme density 



functional with the Skill parameterization of Ref. [29]. Our results for the BCS phase and 
BCS-BEC crossover are consistent with earlier studies: we observe a smooth crossover to 
an asymptotic state corresponding to a mixture of a deuteron Bose condensate and a gas 
of excess neutrons. The transition from BCS to BEC is established according the following 
criteria: (i) The average chemical potential Jx changes its sign from positive to negative 
values, (ii) the coherence length of a Cooper pair becomes comparable to the interparticle 
distance, i.e., £ ~ d ~ p" 1 / 3 as conditions change from £ 3> d to £ <C d. 

The nuclear LOFF phase arises as a result of the energetic advantage of translational 
symmetry breaking by the condensate, in which pairs acquire a non-zero CM momentum Q. 
As illustrated in Fig. [T] at log(p/po) = —0.5, where po = 0.16 fm -3 is the nuclear saturation 
density, the gap in the LOFF phase at non-zero asymmetries and constant temperature 
has its maximum at finite Q, which results in a maximum of the condensation energy 
of the pairs. For large asymmetries the maximum gap occurs for large values of Q. At 
constant asymmetry, a temperature increase shifts the gap maximum and the free-energy 
minimum of the LOFF phase toward small Q, and at sufficiently high temperature and small 
asymmetry the BCS state is favored over the LOFF phase. This behavior is well understood 
in terms of the phase-space overlap of the Fermi surfaces of neutrons and protons, which (at 
finite asymmetry) increases with temperature and the momentum Q of the Cooper pairs. 
Thus, as the temperature increases, we expect a restoration of the BCS phase and of the 
translational symmetry in the superfluid. Obviously, the same restoration occurs when the 
isospin asymmetry is small enough. 

The superfluid phase with phase separation (PS) has the symmetrical BCS phase as one 
of its components. The temperature dependence of this phase is well established within 
BCS theory. The second component, which accommodates the neutron excess, is a normal 
Fermi liquid whose low-temperature thermodynamics is controlled by the excitations in the 
narrow strip of width T '/epn/p around the Fermi surfaces of neutrons and protons. 

The transition to the BEC regime of strongly-coupled neutron-proton pairs, which are 
asymptotically identical with deuterons, occurs at low densities. As already well established, 
in the case of neutron-proton pairing the criteria for the BCS-BEC transition are fulfilled, 
i.e., fx changes sign and the mean distance between the pairs becomes larger than the 
coherence length of the superfluid. 
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FIG. 1: Properties of the nuclear LOFF phase for log(p/po) = —0.5. Left panel: Dependence of 
the pairing gap and free energy of the LOFF phase on the total momentum Q of a Cooper pair 
at T = 1 MeV for asymmetries a = 0.2 (dashed line), 0.3 (dashed-dotted line), and 0.4 (dashed- 
double-dotted line). Right panel: Same dependence as in the left panel, but for fixed asymmetry 
a = 0.2 and temperatures 1 MeV (dashed line), 2 MeV (dotted line), 3 MeV (solid line). 



We now turn to the question of how the BCS-BEC crossover is affected by the existence 
of nuclear LOFF and PS phases at non-zero isospin asymmetries, and conversely how these 
phases evolve in the strongly-coupled regime if the density of the system is decreased. The 
phase diagram of pair-correlated nuclear matter in the density and temperature plane is 
shown in Fig. [2] for several isospin asymmetries. Four different phases of matter are present 
in the diagram, (i) The unpaired phase is always the ground state of matter at sufficiently 
high temperatures T > T c0 , where T c0 (p) is the critical temperature of the superfluid phase 
transition at a = 0. (ii) The LOFF phase is the ground state in a narrow temperature- 
density strip at low low temperatures and high densities, (iii) The PS phase appears at 
low temperatures and low densities, while the isospin-asymmetric BCS phase is the ground 
state for all densities at intermediate temperatures. In the extreme low-density and strong- 
coupling regime the BCS superfluid phases have two counterparts: the BCS phase evolves 
into the BEC phase of deuterons, whereas the PS-BCS phase evolves into the PS-BEC phase, 
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FIG. 2: Phase diagram of dilute nuclear matter in the temperature-density plane for several isospin 
asymmetries a. Included are four phases: unpaired phase, BCS (BEC) phase, LOFF phase, and 
PS-BCS (PS-BEC) phase. For each asymmetry there are two tri-critical points, one of which is 
always a Lifshitz point. For special values of asymmetry these two points degenerate into single 
tetra-critical point at log(p/ po) = —0.22 and T = 2.85 MeV (shown by a square dot) for 04 = 0.255. 
The LOFF phase disappears at the point log(p/po) = —0.65 and T = (shown by a triangle) for 
a = 0.62. 

in which the superfluid fraction of matter is a BEC of deuterons. The superfmid-unpaired 
phase transitions and the phase transitions between the superfluid phases are of second order 
(thin solid lines in Fig. [2]), with the exception of the PS-BCS to LOFF transition, which is 
of first order (thick solid lines in Fig. [2]). The BCS-BEC transition and the PS-BCS to PS- 
BEC transition are smooth crossovers. At non-zero isospin asymmetry the phase diagram 
features two tri-critical points where the simpler pairwise phase coexistence terminates and 
three different phases coexist. (We do not include the points associated with crossovers 
from strong to weak coupling in the class of critical points, since these transitions involve 
essentially the same phase, no symmetry being broken). 

The topology of the phase diagram and the location of the tri-critical points depends on 



8 



the value of asymmetry parameter. For a < the low-density critical point corresponds 
to coexistence of BCS, PS, and LOFF phases, whereas the high-density critical point corre- 
sponds to coexistence of LOFF, BCS, and unpaired phases and is thus a Lifshitz point [30 ]. 
For a > CK4 the topology of the phase diagram changes: the low-density tri-critical point 
contains BCS, PS, and unpaired phases, whereas the high-density tri-critical Lifshitz point 
contains the LOFF-PS-unpaired triple of phases. Clearly, the point with log(p/po) — —0.22, 
T = 2.85 MeV, and = 0.255 is the special case of a tetra-critical point, where all four 
phases, i.e., BCS, PS, LOFF, and unpaired phases all coexist. 

The extreme low-density region of the phase diagram features two crossovers. At inter- 
mediate temperatures we recover the well-known BCS-BEC crossover, where the neutron- 
proton BCS condensate transforms smoothly into a BEC gas of deuterons with some excess 
of neutrons. The new ingredient of our phase diagram is the second crossover at low tem- 
peratures, where the heterogeneous superfiuid phase is replaced by a heterogeneous mixture 
of a phase containing a deuteron condensate and a phase containing neutron-rich unpaired 
nuclear matter. 

Conclusions. Our comprehensive study of pairing in dilute nuclear matter shows that 
at non-zero isospin asymmetry and at low temperatures, pair-correlated nuclear matter 
has a rich phase structure with multiple superfiuid phases and critical points. The phase 
diagram contains three superfiuid phases (BCS, LOFF, and PS) and the strong-coupling 
counterparts (Bose condensates) of two of them. We have called attention to the existence 
of a heterogeneous, phase-separated phase in the context of isosinglet pairing in nuclear 
matter, demonstrating that this phase can occupy a large fraction of the low-temperature 
and low-density region of the phase diagram. We find that the nuclear-matter phase diagram 
contains two tri-critical points where three of the four phases studied coexist; furthermore, 
at a certain value of isospin asymmetry these points degenerate into a single tetra-critical 
point. These diverse features of superfiuid, asymmetrical nuclear matter promise significant 
ramifications for the astrophysics of supernova and (hot) compact stars and therefore warrant 
examination in further detail. 
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